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ABSTRACT 

It is shown that any convolution operator in the time domain can be rep- 
resented exactly as a multiplication operator in the time-scale (wavelet) 
domain. The Mellin transform gives a one-to-one correspondence be- 
tween frequency filters (system functions) and scale filters (multiplica- 
tion operators in the scale domain), subject to the convergence of the 
defining integrals. The usual wavelet reconstruction theorem is a special 
case. Applications to the denoising of random signals are proposed. It is 
argued that the present method is more suitable for removing the effects 
of atmospheric turbulence than the conventional procedures because it 
is ideally suited for resolving spectral power laws. 
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1. Scale Filtering in the Wavelet Domain 

Time-scale analysis of signals is similar to their time-frequency analysis, but with the 
frequency replaced by a scale parameter [1]. Instead of having to choose a basic window 
in time, one must begin with a basic wavelet il^{t). We define the wavelet family of i{j{t) 
as the two-parameter family of functions 

V'cr,T(^) = i^icrt — t) with a, r real and a ^ 0. (1) 

Note that this differs from the usual convention tps^r'it) = \s\~^^'^ip{{t — t')/s). As will 
be seen below, (1) is equivalent but formally simpler. Whereas the usual scale parameter 
s can be interpreted as time scale {\s\ 3> 1 means coarse scale and < |s| <C 1 means 
fine scale) , our parameter a — 1/ s can be interpreted as frequency scale: Large a isolates 
high frequencies, while small cr isolates low frequencies. Thus we may identify a with the 
frequency band passed by the wavelet ifj^^rit), and (a, r) can be regarded as a pair of time- 
frequency parameters, with the understanding that now "frequency" means "frequency 
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scale." That is, instead of shifting {f —>■ f + fo), our frequencies now multiply (a — > aao). 
The absence of the factor = |o"|^/^ will be seen to be unimportant and, moreover, 

to lead to simpler reconstruction formulas. Our parameter r also differs from the above 
r' by T = ar' . That is, r is a time shift in the scaled time. (First scale, then shift rather 
than the other way around.) The wavelet transform of a signal x(t) with respect to the 
family is defined as the inner product 

/oo 
dtij{at-Trx{t). (2) 
-oo 

Inserting the Fourier expansion of 'i/j{at — r) gives 

/oo nOO 
dt / df e-'^'f(^'-^Hifrxit) 
-oo J — oo 

/oo 
df e2-^^-V'(/)*x(/^) (3) 
-oo 

where lipif) denotes the Fourier transform of 'ip{t) and J^~^ is the inverse Fourier transform 
operator. Thus by Plancherel's theorem. 



/oo /"OO 
dT\x{a,T)f= / df\i^{f)f\x{fa)\ 
-OO J — oo 



oo 

oo 



1 r°° 

- / df\^{f/a)\'\x{f)\' 

fI J-oo 



(4) 



Suppose we integrate both sides over a with an arbitrary positive weighting function w{a). 
(For reasons to become clear below, we need all scales a 0. For some purposes, cr > 
will do; this can be arranged simply by choosing w{a) = for u < 0.) Thus 



/oo POO poo 1 poo 

daw{a) / dr |x(^,t)|2= / —w{a) / df mf /a)\' \x{f)\ 
-oo J — oo J —oo \^\ J — oo 



oo 1^ I J — oo 

oo /"OO 



/oo pOO J 

df\m? ^w{a)m/a)\\ 
-oo ./ — oo I I 



(5) 



where the integral over cr is understood to exclude (7 = (integrate over < £ < |(t| < oo, 
then take the limit £ — > 0). Let 

*(/) = m)? 

denote the spectral density of '^{t), and define 

/oo 7 
— w{a)^>{f/a)^{w.^>){f) 
-oo l-l 

^{a)w{f/a)^{^.w){f\ 



oo 
oo 
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where the second hne is obtained from the first by the substitution a — > //cr, assuming 
that f ^ 0. (The "DC component" / = gets special treatment throughout wavelet 
theory, as does (7 = 0, and for similar reasons.) Then (5) becomes 

j I dadrwia) \x{c7,t)\' = df W{f) |x(/)p. (7) 

The left-hand side in (7) is a scale-weighted norm ("energy") of the signal, whereas the 
right-hand side is a frequency-weighted norm with weight function W{f). (For example, 
Sobolev norms are defined with W{f) = (1 + f^Y for some power s.) Then (7) is a kind 
of "Plancherel theorem" equating the two weighted norms. Like the Plancherel theorem, 
it can be polarized to give the analog of Parseval's relation for two signals 0(t) and x(^)- 

JJdadT 0(a, T)*«;(a) xia, r) = J^ df 0(/)*W-(/) x(/)- (8) 

This shows that the convolution operator W with system function (or "symbol") W{f) 
can be expressed in the wavelet domain as 

(Wxm = J^ df e^-'f'Wif)xif) = JJ dadr iP^^t) wia) xia^r). (9) 

Thus, to represent W in the waveJet domain, we need merely to replace x(/) with x{a, r), 
W{f) with w{a), and the "basis vectors" e^"^*-^* of the Fourier transform with the "basis 
vectors" V'a,r(^) of the wavelet transform! 

In deriving (9) we assumed that w{a) (and therefore, by (6), also W{f)) is nonnegative. 
This assumption was made for technical reasons, since then all the integrands occurring in 
Equations (5)-(8) are nonnegative and the integrals therefore either converge absolutely 
or diverge to infinity, making it unnecessary to worry about conditional convergence. We 
now relax this assuption and allow w{a) (hence also W{f)) to be complex-valued, with the 
understanding that the integrals may no longer converge absolutely. Then (9) shows that 
operations normally performed in the frequency domain can now also be performed in the 
time-scale domain, provided we can find a weighting function w{a) which gives W{f) by 
(6). So the first problem that must be addressed is: Given W{f), solve (6) for w{a). We 
now show that this problem has a unique solution, subject to the admissibility condition 
that its defining integrals converge. 

The key is to note that Equation (6) defines W{f) as a scaling convolution of w{a) with 
^{f), in the sense that the usual difference variable / — cr is replaced by the quotient //cr 
and the translation-invariant (Lebesgue) measure da is replaced by the scaling-invariant 
measure da/\a\. Just as convolutions are converted into products by the Fourier transform, 
scaling convolutions are converted into products by the Mellin transform [2]. The latter 
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takes a function F{a) of a positive variable cr > into a function F{p) of a complex 
variable p by 

poo 1 

M{F{a)}{p) = —a-^ F{cj) = (10) 
Jo 

and the original function is reconstructed by integrating over a contour C in the complex 
plane of the "power parameter" p going from c — ioo to c + ioo, for an appropriate choice 
of c which depends on F{(t): 

Fia) = 1-Jjpa^ Fip) = M-'{Fip)}ia). (11) 

To solve (6) for w{a) using the Mellin transform, we must separate the positive- and 
negative-frequency components of W{f), since the Mellin transform uses only a > 0. Let 
us assume that 

^(f) = for / < 0, hence *(/) = for / < 0. (12) 

Then the wavelet 'ip{t) is necessarily complex. It is an analytic signal in the sense of Gabor 
[3], extending analytically to the upper half of the complex time plane. (See also [1], 
Section 9.3.) Let U{a) be the unit step function 



U{a) 



1 if a > 
if CT < 0, 



and write 
where 



w{a) = U{a)w{a) + U{-a)w{a) = w+{a) + w_{-a), (13) 



w±{a) = U{a)w{±a) = for a <0. (14) 
By (12), (6) becomes (for / 0) 

W{f)= r —n'^)w{f/a) = W4f) + W_{-f), (15) 



where 



poo 1 

WM)= / — *(^)«;±(//^) = (*««;±)(/) = Ofor/<o. (16) 

Jo cr 

Therefore the positive-frequency component W+(/) of W{f) depends only on the positive- 
scale component w_^_{a) of w{a) and the negative- frequency component W_{—f) depends 
only on the negative-scale component w_{—a). This decoupling was the purpose of the 
assumption (12) above. Applying the Mellin transform to the scaling convolutions (16), 
we obtain 

W±{p) = • w±] {p) = ^{p) w±{p). (17) 



This gives a formal expression for w±{a) in terms of the transform of W±{f): 



w±{<j) = dpa P a>0. 



(18) 



'ip) 

We give three examples. First, let n be an integer and 
Then 

w^{a) = U{a){±ar, (19) 

and (16) gives 

poo 7 

WM) = U{f) {±fr / — *(a)c7-" = U{f) (±/)"*(n). (20) 

Hence 

W{f) = *(n) [Uif) r + Ui-f) n = *(n) r • (21) 

That is, pure integer powers are invariant under the correspondence w{a) ^ W{f) between 
the frequency and scale domains, except for the renormahzation constant ^(n). The 
differential operator 

P(i?) = ^a„i?^ where D = 

n=0 

is represented in the frequency domain as multiplication by 

AT 

W{f) = J2^n {2nifr. (22) 

n=0 

Hence, by (21), it is represented in the wavelet domain as multiplication by 

. ^ an(27ri(T)'" 
^(^) = 2^ , • (23) 

In order for this representation to be defined, we must have 

< *(n) = / — *(c7) (j"" < oo for allnwithon ^ 0. (24) 
Jo cr 

We say that the operator P{D) is admissible if the condition (24) holds. In particular, the 
identity operator P{D) = 1 is represented in frequency by W{f) = 1, hence 

w{a) 
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and (9) becomes 

Xit) = J^^II dadT^,,r{t) x(^T, r). (25) 

This is the usual wavelet reconstruction formula, which inverts the wavelet transform. It 
shows that the identity operator is admissible if and only if 

*(0)^ / — ma)f<oc, (26) 
Jo o- 

which is the usual admissibility condition [1] for the wavelet '4>{t). Indeed, (18) indicates 
the following general admissibility condition for convolution operators: 



In order for the convolution operator with system function W{f) to be 
representable in the wavelet domain of the wavelet ip{t), the Mellin trans- 
forms W±{p) of VF±(/) must be analytic along an appropriate contour 
C parallel to the imaginary axis such that ^'(p) is also analytic with no 
zeros along C and the integral (18) converges. 



This example suggests that wavelet methods may be used instead of Fourier methods to 
solve differential equations. To make this precise, one must examine the admissibility 
condition of the given operator and the range of the wavelet transform, i.e., the set of 
all functions x{a, r) obtained from x(t) as the latter is allowed to range over the space of 
functions of interest. (Not every function F{a, r) of time and scale is the wavelet transform 
of some time signal x(t); F{a,T) must satisfy a consistency condition related to the fact 
that the "basis vectors" 'il)(j,T{t) are redundant [1].) 

Our second example is 

w{a) = \a\P, (27) 
where p is now any power (not necessarily an integer) . Then 

w±((j) = f/((j)(jP, (28) 

and (16) gives 

WM) = U{f)F r —^{a)a-P = Uif)r^ip). (29) 
Jo o- 

Hence 

w{f) = *(p) [u{f) r + ui-f) i-m = m i/r- (30) 

Like the pure integer powers a'^ of a, the pure powers |cr|^ of |cr| are invariant under the 
correspondence w{a) ^ W{f), again with the renormalization factor ^'(p). Note that this 
representation of the filter W{f) = |/|^ (with p = —a) is implicit in [1], Section 9.3. 

Our final example is 

h{a) = J-\-' (31) 
^ ^ ^(0) li if (7 < 0. ^ ^ 
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This gives 



BU) = {.' mi (32) 



which is the system function of the Hilbert transform 

1 poo 7 

{nxm = -p -^xW, (33) 



where P denotes the principal part of the integral. Thus 



mm 



*(0) 



JJa<0 JJa>0 



(34) 



This representation of Ti is implicit in Theorem 3.3 of [1]. 

The combination of Mellin and Fourier transforms has long been used by harmonic 
analy lists in connection singular integral operators, in particular to obtain estimates [4]. 
But although the concrete representation of convolution operators derived here would seem 
to be fundamental in the application of wavelet analysis to signal processing, I was unable 
to find its existence in the literature. A discrete version of this problem in the framework 
of multiresolution analysis was the subject of a recent paper by Beylkin and Torresani [5] . 



2. Applications 

The time-scale representation (9) can be applied, in principle, to any problem involving 
convolution operators. Its usefulness depends on the nature of the system function W{f) 
for the problem at hand. The above examples show that systems whose functions obey 
power laws will have simple representations. In particular, the problem of filtering out 
the effects of atmospheric turbulence (denoising) seems a likely candidate, since it involves 
asymptotic power laws. There is strong empirical evidence that denoising in the scale 
domain is superior to denoising in the frequency domain because it tends to smooth out 
the spectrum without smearing the scales, so that the exponents in the power laws are 
clearly resolved [6]. 

The above time-scale analysis was developed for one dimension (time). It can be 
extended to several dimensions, like space-time, and in that case it may be applied to 
physical wavelet representations. Physical wavelets, which were defined in [1], are local- 
ized acoustic or electromagnetic waves (i.e., solutions of the scalar wave equation or of 
Maxwell's equations) which share many of the properties of the one-dimensional wavelets. 
They are related to one another by geometrical operations such as space-time translations, 
scalings, and Lorentz transformations (boosting to moving coordinate systems), and they 
can be used as "building blocks" to form arbitrary acoustic or electromagnetic waves. 
Such wavelets have recently found natural applications in sonar and radar [7] . With them. 
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the problem of denoising atmospheric turbulence can be approached from first principles, 
since the physical wavelets already fully describe the deterministic space-time behavior 
(dynamics) of the mean process. Progress on this will be reported elsewhere. 
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